First principles study of electronic transport through a Cu(lll)|graphene junction 
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We report first principles investigations of the nonequilibrium transport properties of a 
Cu(lll)|graphene interface. The Cu(lll) electrode is found to induce a transmission minimum 
(TM) located — 0.68eV below the Fermi level, a feature originating from the Cu-induced charge 
transfer resulting in n-type doped graphene with the Dirac point coinciding with the TM. An applied 
bias voltage shifts the n-graphene TM relative to the pure graphene TM and leads to a distinctive 
peak in the differential conductance indicating the doping level, a characteristic not observed in 
pure graphene. 
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Since the discovery of graphene (Gr), much research, 
fueled by the advancement of scientific and technological 
progress at the nanoscale, has focused on the exceptional 
intrinsic electronic character of this two-dimensional ma- 
terial pQ. Within the context of nanoelectronics, utilizing 
Gr as an element in a device, e.g., electrode or channel, 
most likely requires interfacing it with other materials, 
a process known to potentially override /limit the unique 
properties of the Gr. In this view, recent experimental 
[2] and theoretical [3HS] studies have focused on the con- 
tact properties of Gr-based devices and how the presence 
of other materials can influence the global response of a 
device; a principle used to demonstrate a high-speed Gr 
photodetector [7]. Thus in regard to transport, a physi- 
cal understanding of the Gr- metal contact, including the 
materials specific bonding, is paramount. 

With the advent of large-scale epitaxially grown Gr on 
Cu [8], relatively straightforward fabrication of macro- 
scopic Gr samples is possible, hence facilitating scien- 
tific and technological purposes. This advancement also 
highlights the current importance of Cu-Gr interaction. 
Whether etching part of the Cu substrate to realize a de- 
vice or transferring Gr to another substrate, some degree 
of Cu may remain in contact with the Gr. Thus a clear 
understanding of the influence of Cu on the electronic 
and transport characteristics of Gr is advantageous. 

In this paper, we present a nonequilibrium ab initio 
study of the transport properties of a Cu-Gr interface 
(as shown in FigjTJa), hereafter referred to as Cu|Gr). 
Our findings show that a Cu electrode induces a trans- 
mission minimum (TM) located at — 0.68eV below the 
Fermi level (Ep), denoted as TM1, in addition to the 
TM of pure Gr at E F (TM2), both due to a vanishing 
density of states (DOS) at the Dirac point. TM1 orig- 
inates from a weak hybridization between Gr and Cu, 
which (i) preserves Gr's linear dispersion bands around 
the K point and (ii) induces a net charge transfer from 



*Electronic address: maassenj ©physics. mcgill.ca 
t Electronic address: wji@ruc.edu.cn 



Cu to Gr resulting in n-type doped Gr (n-Gr) [3|. Apply- 
ing a bias shifts the position of the TMs relative to each 
other, leading to a peak in the differential conductance 
(dl/dV) that is otherwise not observed in pure Gr. 

The atomic coordinates of the Cu|Gr interface were 
optimized using VASP [9], a density functional theory 
(DFT) package. The nonequilibrium transport proper- 
ties were calculated using MatDCal [10 , which com- 
bines nonequilibrium Green's functions (NEGF) with 
DFT to model open systems in a two-probe geometry. 
For MatDCal, we use an optimized double-^ polar- 
ized atomic basis set, the local density approximation for 
exchange-correlation energy [11], and norm-conserving 
non-local pseudopotentials [12] to account for the nu- 
clear and core potentials. Converged /c-meshes, including 
the important T, K and M points, were well tested and 
adopted in all calculations. 

As illustrated in Figjlja) the left and right leads, ex- 
tending to ±oo, are comprised of Cu-sandwiched Gr, sim- 
ulating the situation where a metallic contact is deposited 
on Gr, and pure Gr, respectively. The system was chosen 
to replicate the interface between a source/drain region 
of a Gr device (Cu-coated Gr) and the pure Gr channel. 
Periodic boundary conditions are imposed in the plane 
perpendicular to transport. Figure[TJb) shows the atomic 
structure of the left lead of the Cu(lll)|Gr interface. The 
most stable configuration of Gr on Cu(lll) corresponds 
to one sub-lattice C atom situated directly above a Cu 
atom and the other C atom located at the hollow site 
of the Cu(lll) surface [3 j. The supercell of the left lead 
for total energy and electronic structure calculations is 
shown with the dotted box in FigjTJb), where the images 
of the Gr sheets in neighboring supercells are separated 
by seven Cu(lll) layers. By simultaneously relaxing the 
atomic coordinates to forces < 0.01 eV/A and optimizing 
the supercell box size by total energy minimization, the 
Gr-Cu(lll) interlayer distance is found to be 2.95 A. 

The band structure of the left lead is shown in Fig[TJc). 
The blue lines represent the hybrid states of the Cu-Gr 
contact, while the black circles show the p z character of 
the C atoms corresponding to the Gr 7r-states. From this 
projection, we observe that the Gr bands conserve a lin- 
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FIG. 1: (Color online) (a) Atomic structure of the Cu(lll)|Gr 
interface. The black dotted lines delimit the scattering region 
from the left lead (Cu-contacted Gr) and right lead (pristine 
Gr) each extending to =boo. (b) Side- and top- view of the left 
lead. The optimized distance between the Gr sheet and the 
Cu(lll) surface is 2.95 A. (c) Band structure of the left lead. 
The states originating from the hybrid Cu(lll)|Gr system are 
plotted as solid blue lines, while the C(p z ) character of the 
bands are superimposed as black circles. The dotted black 
lines, shown in (b), indicate the supercell box used for the 
electronic structure calculation. 



ear dispersion near the K point, as previously shown for 
Gr in contact with one Cu(lll) surface [3]. The weak 
hybridization between Gr and Cu merely shifts the con- 
ical Gr states down — 0.68eV relative to Ep resulting in 
n-type doped Gr. This energy shift in the Gr bands 
is somewhat larger than the value obtained for a single 
Cu(lll) surface with the same interlayer distance [3]. 

Figure [2] shows the transmission coefficient (T) versus 
energy (E) as a function of Gr orientation at equilibrium 
(y = 0). A Gr sheet oriented such that transport oc- 
curs along the zigzag or armchair direction is denoted as 
Gr (zigzag) or Gr (armchair), respectively. The solid thick 
blue (thin red) line corresponds to the Cu|Gr (zigzag) 
(Cu|Gr(armchair)) interface. A TM (i.e., TM2) is ob- 
served at Ep (Ep is shifted to zero) , due to the vanishing 
density of states (DOS) in the pure Gr right lead. Inter- 
estingly, a second TM (i.e., TM1) is located at -0.68 eV. 
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FIG. 2: (Color online) Transmission versus energy at equi- 
librium (V — V) as a function of Gr orientation. The 
solid thick blue (thin red) line represents a Cu(lll)|Gr in- 
terface where the Gr is oriented such that transport occurs 
along the zigzag (armchair) direction. The dashed thick blue 
(thin red) line shows the transmission for pure zigzag-oriented 
(armchair-oriented) Gr. The Fermi level is set to zero. 



This feature also originates from the near-zero DOS at 
the neutrality point of the Cu-induced n-Gr in the left 
lead (see Fig[]Jc)) [6]. This indicates that the many avail- 
able metallic bands near the TMs, shown in the band 
structure, play a negligible role in the electron flow of 
the junction. To study the influence of the Cu electrode 
on the transport properties, we calculate T for pure Gr 
(dashed lines). We find that T near Ep is unaffected by 
the presence of the Cu electrode, while the same is also 
true near the TM1. This confirms that a Cu(lll)|Gr in- 
terface can be effectively viewed as a n-Gr|Gr junction. 
Finally, we consider the effect of Gr orientation on T 
and conclude that both curves are related by a constant 
scaling factor. The ratio of T z i gzag /T armc h a ir is found to 
be equal to the ratio Warmchair/Wzi gzag , where W is the 
width of the supercell used for transport. Thus, T (and 
hence current) per unit width results in identical values 
for both zigzag and armchair orientations. 

In Fig{3^a), we present T versus E as a function of 
applied bias (V) for Cu|Gr (zigzag). A bias V varies the 
chemical potentials of the left and right leads as hl = Ep 
and jiR = Ep + |e|V, where e is the electron charge. 
From Fig(3^a), we see that the TM1 is pinned at —0.68 eV 
(shown with the dashed vertical line), whereas the TM2 
moves by an amount \e\V (indicated by the arrows). This 
is easily understood because our bias is applied to the 
right Gr lead which shifts the Dirac point there. By 
applying a bias one can shift the relative positions of 
the TM features. From the T calculated at finite V, we 
obtain the current per unit width (I) for Cu|Gr (zigzag) 
and pure Gr (zigzag) shown in Figj3[b). I is calculated 



3 
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FIG. 3: (Color online) (a) Transmission versus energy as 
a function of applied bias V for the Cu(lll)|Gr (zigzag) 
interface. The transmission curves corresponding to 
V = +0.8 V, +0.4 V, V, -0.4 V, -0.8 V, are offset by 
0, 0.2, 0.4, 0.6, 0.8 respectively, with the zero plotted as a 
dashed horizontal line. The arrows (black dashed vertical 
line) indicate the position of the Dirac point of the pure Gr 
right lead (n-type Gr left lead). Current / (b) and differ- 
ential conductance dl/dV (c) per unit width versus V for 
Cu(lll)|Gr(zigzag) and pure Gr (zigzag). 
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where h is Planck's constant, W the unit cell width and 
/l,r the Fermi-Dirac distribution functions. The I-V 
curves are smooth non-linear functions of V. To em- 
phasize the differences between the Cu|Gr and pure Gr, 
we calculate the differential conductance per unit width 
G = dl/dV, presented in Figj3^c). Pure Gr results in a 
symmetric monotomic G, whereas Cu|Gr shows a peak 
near —0.63 V where both Dirac points from the Gr and 
n-Gr coincide in energy. For V greater than —0.63 V, the 
TM1 plays no role in the values of / and G, since it is 
out of range for the integration in Eq.Q (from E F + \e\ V 
to Ep, when V < and temperature is zero). When 
V < —0.63 V, the TM1 enters the energy integration 
range for / and results in a rapid decrease in G. The 
maximum derivative of G provides the location of the 
conical point of the n-Gr. This prediction should be ex- 
perimentally detectable. 

In summary, we performed first principles transport 
studies of the Cu(lll)|Gr interface. A Cu(lll) elec- 
trode on Gr induces a second conductance minimum near 
— 0.68eV relative to Ep, originating from the vanishing 
density of states of the Cu-induced n-doped Gr. We find a 
Cu|Gr interface can be effectively modelled as a n-Gr|Gr 
junction due to a relatively weak interaction between Cu 
and Gr. An applied bias shifts the positions of the con- 
ductance minimums relative to each other, leading to a 
peak in the differential conductance (indicating the dop- 
ing level of the n-Gr), a distinctive feature not observed 
in pure Gr. These results provide a general picture of 
the influence of weakly interacting metallic electrodes, 
including Al, Ag, Au and Pt [13], on transport in Gr. 
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